Dopant-modulated pair interaction in cuprate superconductors 
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Comparison of recent experimental STM data with single-impurity and many-impurity 
Bogoliubov-de Gennes calculations strongly suggests that random out-of-plane dopant atoms in 
cuprates modulate the pair interaction locally. This type of disorder is crucial to understanding 
the nanoscale electronic structure inhomogeneity observed in BSCCO-2212, and can reproduce ob- 
served correlations between the positions of impurity atoms and various aspects of the local density 
of states such as the gap magnitude and the height of the coherence peaks. Our results imply that 
each dopant atom modulates the pair interaction on a length scale of order one lattice constant. 

PACS numbers: 74.72.-h,74.25.Jb, 74.20.Fg 
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The discovery of nanoscale inhomogeneity in the 
cuprates has recently generated intense interest. In par- 
ticular, the spectral gap in the local density of states 
(LDOS), as observed by scanning tunneling microscopy 
(STM) QI1HI3 in Bi 2 Sr 2 CaCu 2 08 +a; (BSCCO), varies 
by a factor of two over distances of 20-30A. This un- 
usual behavior may help reveal how the cuprates evolve 
from the Mott insulating state at half-filling to the su- 
perconducting state at finite doping. The hole concen- 
tration in the Cu02-planes of BSCCO is proportional to 
the number of out-of plane dopant atoms, which also in- 
troduce disorder. This has led to the proposition that 
poorly screened electrostatic potentials of the dopant 
atoms generate a variation in the local doping concentra- 
tion and thus give rise to the gap modulations observed 
in STM [H, E3, ■ Poor screening has also been argued to 
result in enhanced forward scattering |8| , which appears 
to be compatible with photoemission jij [ljj and trans- 
port measurements in the superconducting state of 
BSCCO. An alternate perspective is explored in several 
works which associate inhomogeneous electronic struc- 
ture with a compet ing order parameter, such as anti- 
ferromagnetism [l2l fl3L llij . Only very recently has it 
been possible to measure correlations between the inho- 
mogeneities observed in STM and positions of dopant 
atoms , thus providing a clue to the relation between 
disorder and doping in this compound, as well as a means 
to examine the above proposals. 

In this Letter, we assume that the electronic inhomo- 
geneity observed by STM, at least in the optimally to 
overdoped samples, can be understood within the frame- 
work of BCS theory in the presence of disorder. We show, 
however, that the conventional modeling of disorder as 
a set of random potential scatterers fails to reproduce 
the most prominent features of the STM experiments 
described above: (i) the subgap spectra are spatially ex- 
tremely homogeneous Q, unless they are taken in the 
immediate vicinity of a defect in the Cu02 plane, (ii) 
the coherence peaks in regions with larger gap tend to 
be much broader and reduced in height, (iii) the "co- 
herence peak" positions are symmetric about zero bias, 



(iv) the dopants are found to correlate positively with 
large gap regions [lij : and (v) charge modulations are 
small We propose that the dopant atoms mod- 

ulate the local pair potential, i.e. the local attractive 
coupling g between electrons is spatially dependent. In 
conventional superconductors, such effects are difficult to 
observe because atomic-scale modulations in g produce 
LDOS modulations only on the scale of the coherence 
length £o- In the cuprates, however, the situation is dif- 
ferent due to the short coherence length. Wc demonstrate 
that a model in which dopant atoms modulate the pair 
interaction gives excellent agreement with respect to the 
above mentioned key characteristics of the STM data. A 
modulated pair interaction could arise from local lattice 
distortions surrounding the dopant atoms modifying the 
electron-phonon coupling or superexchange interaction in 
their vicinity. 

Model. We consider the following mean-field Hamilto- 
nian for a singlet d-wave superconductor 



(1) 



where = —2t(cosk x + cosk y ) — At 1 cos k x cos k y — /1 
and denotes summation over neighboring lattice 

sites i and j. In the remainder of the paper we will set 
t' /t = —0.3 and adjust fi to model the Fermi surface 
of BSCCO near optimal doping (for the homogeneous 
system, /Lt/t=-1.0). In order to account for disorder in 
the out-of-plane dopants, which are separated from the 
Cu02 plane by a distance z, we include an impurity po- 
tential modeled by Vi = Vq exp(— rj/A)/r,, where r, is the 
distance from a dopant atom to the lattice site i in the 
plane. Distances are measured in units of V2a, where 
a is the Cu-Cu distance. The nearest-neighbor d-wave 
order parameter = gij(ci]Cj[ — CjiCii) is determined 
self-consistently using with g^ = g + (Vi + Vj)/2, In 
traditional BCS theory, <?y = g is spatially uniform, and 
Aij is only modulated in the vicinity of potential scatter- 
ers [la, LLZj . We will argue that this approach is unable to 
reproduce observations (i)-(v) outlined above, and that 
gij is strongly modified near the dopant atoms. 
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Smooth potential. If the potential caused by the out- 
of-planc dopant atoms were very smooth on the scale 
of £ , the local properties of the inhomogeneous system 
would be determined by the local value of the disor- 
der potential and the local value of the pairing inter- 
action. Therefore one would expect an LDOS which is 
locally similar to a clean superconductor with renormal- 
ized chemical potential [m — Vi in the case of a smooth 
potential Vi, or with renormalized bond order parame- 
ter Aij + SAij for a smooth off-diagonal (OD) potential. 
In the case of a conventional diagonal potential, a gap 
size modulation will be induced because the gap is a rel- 
atively sensitive function of the local chemical potential, 
see Fig. QJa). On the other hand, modulations of this 
type will inevitably have coherence peak weight-position 
correlations opposite to experiment, since large gap val- 
ues in the homogeneous system imply (within BCS the- 
ory) that spectral weight removed from low energies is 
transferred into the coherence peaks (Fig. [U[b)). This 
effect is further enhanced by the presence of a van-Hove 
singularity at lo v h = s f{U l + v) 2 + (4A ) 2 in the tight- 
binding model which contributes additional weight to the 
coherence peaks, in particular for fi/t = —1.2 where it co- 
incides with the gap edge. Here An. is the bond order pa- 
rameter in the homogeneous system. A similar although 
less pronounced positive correlation between coherence 
peak weight and position arises also for the smooth OD 
case. Note that throughout this work we neglect inelastic 
scattering that would broaden the tunneling conductance 
peaks at large bias but would not change their weight, 
thus leaving our conclusions unaffected. 

Single-impurity scattering. Since a smooth disorder 
potential cannot reproduce the experimentally observed 
relation between the weight of the coherence peak and 
the gap magnitude, we now address the opposite limit, 
i.e. a very spiky potential caused by a dopant poten- 
tial with short range on the scale of £n. Some insight 
into this situation can be obtained by analyzing single- 
impurity scattering processes, which should be dominant 
for sufficiently short ranged and weak scattering potcn- 
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FIG. 1: (a) Order parameter Ao for t' /t = —0.3 and a con- 
stant nearest neighbor attraction of g/t = 1.16 as a function of 
chemical potential fj,. (b) Local density of states for different 
a using the order parameter displayed in (a). 





-i 1 


1 1 1 1 1 V 


-\ 


1 i 1 i 1 i- 


-j 1 


1 1 1 1 1 I 


:y site 


(a) 


— v=o 

V/t=0.4 - 

--■ V/t=-0.4 


(b) 


— bulk 

4 bonds 

--■ extended 


(c) 


! — bulk 

J 4 bonds 

\ — extended - 

j, ( 














!' * 

l h 


impui 


_ 










i ii 

j! ! 












i 'i 














i 'i 


Ldos 


i , 






X>*r i , i 


1 , 1 


■ 1 1 ' Im 

Vi i'/i/A 
J i.' J -. V 

v — 

A/ 1 . 1 



-0.8 -0.4 0.4 0.8 -0.4 0.4 -0.8 -0.4 0.4 0.8 



CO/t CO/t CO/t 

FIG. 2: On-site LDOS for different single-impurity models 
with t'/t = -0.3, u/t = -1 and A /t = 0.1. (a) Weak point- 
like potential scatterer. (b) Dotted line: attractive "pointlike" 
OD scatterer with 5 A = Ao on the four bonds surrounding 
the impurity site. Dashed line: extended attractive OD scat- 
terer with A = z = 1. (c) Same as (b), with SA — —Ao. 



tials, where interference effects are negligible. For sim- 
plicity we assume constant order parameter in the fol- 
lowing T-matrix analysis and postpone the fully self- 
consistent treatment to the many impurity case. 

Although solving the T-matrix equation for the poten- 
tial scatterer Vi — VSio is straightforward , relatively 
little attention has been paid to the weak to intermedi- 
ate strength impurity case V < t because it does not lead 
to well-defined resonant states inside the gap. Fig. Et a ) 
shows the LDOS at the impurity site for a weak pointlike 
potential scatterer. The positions of the coherence peaks 
are hardly shifted at all, and while the spectral weight 
of the coherence peaks is modified, this occurs in a dis- 
tinctly particle-hole asymmetric fashion. This is in strik- 
ing contrast to the STM spectra, where inhomogeneous 
but particle-hole symmetric coherence peak modulations 
are observed. In addition, there is no distinct feature 
in experiment corresponding to the van Hove features 
present, as e.g. in Figs-^b) and Eta). 

These shortcomings of the conventional potential scat- 
tering model can be overcome by considering OD scat- 
tering instead. For the sake of clarity, in this para- 
graph we neglect the diagonal component of the poten- 
tial. Fig. Etb) and Fig. EJc) show the LDOS at the im- 
purity site for a "pointlikc" OD scatterer with d-wave 
symmetry on the four bonds emanating from site i = 0, 
£An,±£ = — SAo t ±y = SA, and a more extended OD scat- 
terer with SAij = ±SA(Vi + Vj)/(2V ), where Vi is de- 
fined below Eq. JTJ and the negative sign applies to bonds 
oriented along the y-direction. Scattering at an order pa- 
rameter enhancement (see Fig.EJb)) strongly suppresses 
the coherence peaks for large values of SA or more ex- 
tended OD scatterers. For scattering by a local order 
parameter suppression (see Fig. Etc)), exactly the oppo- 
site happens: an Andreev resonance forms just below the 
gap edge, similar to the case where the order parameter is 
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suppressed near surfaces [19|, |20j . For large negative val- 
ues of (5A, or more extended OD scatterers, the Andrecv 
resonance moves to smaller energies, and its peak height 
increases. It draws most of its spectral weight from the 
van Hove singularity at (7r,0), which is close to the part 
of the Fermi-surface with the largest <i-wave gap, i.e., the 
part which is most affected by order parameter modu- 
lations. Although this indicates that the weight of the 
resonance depends on band structure, we find that the 
phenomenon is very robust over a wide range of t' and /i. 

Many -impurity results. In order to strengthen our 
argument for dopant-modulated pairing interaction we 
now address the effect of self-consistency and interfer- 
ence between many impurities. To this end, we solve 
self-consistently the Bogoliubov-de Gennes (BdG) equa- 
tions resulting from Eq.0J, on a 80 x 80 lattice rotated 
by 7r/4 compared to the Cu-0 bond direction (as in ex- 
perimental STM maps), i.e., our system contains 2 x 80 2 
lattice sites in total. Assuming the dopants are intersti- 
tial oxygens, each one most likely contributes two holes to 
the Cu02 plane. We therefore consider a dopant concen- 
tration of 7.5% for optimal doping which are distributed 
randomly in the reservoir layer separated from the CuC>2 
plane by a distance z. 

In the limit of a smooth potential (Fig. EJa)), the 
many-impurity results agree well with the local \i picture 
discussed above. The correlation between the dopant 
positions and the gap amplitude depends strongly on 
the size of the potential due to the non-monotonic de- 
pendence of A on the local as shown in Fig. ^a). 
The spatial variation of the gap, however, is not rapid 
enough to reproduce the grainy gap maps seen experi- 
mentally with gap "patches" of typical size 20-30A 0; 
one is therefore forced to consider "spikier" potentials 
(Fig. Bib)). In the weak limit V < t, one recovers 
the results of the single-impurity case, i.e., the coher- 
ence peaks are modulated in a particle-hole asymmetric 
way. For the stronger spiky potentials required to repro- 
duce the magnitude of the gap modulations observed in 
STM, subgap states start to form in contradiction with 
experiment (sec |3] (b)). Further discrepancies between 
Figs. E2a) and &b) and the experimental spectra are: 
i) the LDOS clearly does not exhibit the inverse rela- 
tion between gap size and coherence peak height; ii) the 
spectra are quite particle-hole asymmetric (see Fig. |^b) 
and high energy regions in Fig.[3fa)); and iii) the sizable 
potential required to induce gap modulations inevitably 
leads to large (0(50%— 100%)) local charge modulations. 
The latter point puts strong constraints on any potential 
scattering model, since the primary role of the impurity 
potential is to couple to the density. 

A typical LDOS line scan for a many OD impurity cal- 
culation is shown in Fig.^Jc). Note that, by construction, 
this model has homogeneous low-energy LDOS as well as 
strong correlations between the dopant positions and the 
local gap values. As in the single impurity case, the linc- 
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FIG. 3: LDOS from self-consistent solution of BdG equations, 
along a straight line for (a) conventional potential with z = 2, 
A = 2, Vo = 1.5t; (b) same as (a), but with z = 0.57, A = 0.5; 
(c) OD potential with Vi as in (b); and (d), combination of 
OD potential shown in (c) with conventional potential as in 
(b) but with Vo = 0.6t. Conventional (OD) potentials are 
depicted to right of each panel as a thin (thick) line. 

shape of the LDOS near the gap edge is determined pri- 
marily by Andreev scattering. Because the LDOS near 
the gap edge is reminiscent of a coherence peak, we will 
simply adopt this terminology, as used in experiment. 
The dopant atoms inevitably give rise to a conventional 
potential as well, however; we therefore show in Fig. Et<d) 
that the qualitative features of OD scattering depicted 
in Fig.^c) survive this scattering. Comparing Figs.^a- 
d), it is evident that the OD LDOS spectra are far more 
particle-hole symmetric than those with potential disor- 
der, and display the inverse relation between gap mag- 
nitude and coherence peak height, as expected from the 
single-impurity discussion (see Fig- EJ) - In Fig. 0] we show 
the associated gap map (a), the coherence peak height 
map (b), and the charge modulation map (c) for pa- 
rameters corresponding to Fig. Ofc). Fig. QJd) displays 
the correlation functions between the gap map and the 
dopants, and the gap map and the peak height map |2lj . 
The local pairing modulation shown in Fig. ^reproduces 
qualitatively the correct negative correlation between the 
gap amplitude and the coherence peak height, the posi- 
tive correlation between dopant atom locations and large 
gap values, and the relatively small charge modulations 
observed in experiment [l5j . In addition, the spectra ex- 
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FIG. 4: (Color online.) Many OD impurity model for param- 
eters of Fig.EJc): (a) 2D real -space map of the local coherence 
peak position (gap) in units of t; (b) coherence peak height 
(note the inverse color scale with respect to (a)); (c) total 
charge (note the small scale) ; and (d) the correlation function 
between the gap map and the dopant atoms (solid), and the 
gap map and the peak height map (dashed). Lines without 
(with) symbols correspond to the parameters used in Fig.^Jc) 
(Fig.Hd)). 



hibit the same remarkable particle-hole symmetric modu- 
lations of the coherence peaks observed in experiment Q . 
This symmetry should manifest itself in Fourier trans- 
form quasiparticle interference patterns as well. 

In the OD scattering model, short-distance correlations 
between the dopant atoms and the gap size are nearly 
perfect, as seen in Fig. 2fd); indeed, they are consid- 
erably stronger than reported in experiment |l5l |. This 
might be due to the difficulty of identifying all dopant 
positions experimentally, to the presence of additional 
cation disorder in BSCCO [2^. or to the finite experi- 
mental resolution of the dopant resonances. The dopant- 
gap correlations are quite robust against inclusion of a 
conventional scattering component (the two solid curves 
in Fig. Etd) coincide), but the gap-peak height correla- 
tions are rapidly suppressed, as seen in Fig.^d). 

A natural question is the extent to which these cor- 
relations are robust against different choices of param- 
eters. We find that the local spectral properties in the 
spiky regime of the OD model are insensitive to param- 
eters, provided the amplitude of the gap modulation SA 
is comparable to or larger than the splitting of the van 
Hove and coherence peaks in the pure system. In that 
case the weight of the van Hove peak is absorbed into the 
coherence peak (Fig. Etc)). 

While we assert the primacy of the OD channel of scat- 
tering for the modulation of the states near the antinodc, 



we emphasize that nodal quasiparticlcs are very weakly 
scattered by this potential, and so microwave and ther- 
mal transport are probabl y on ly minimally influenced by 
the effects discussed here This further implies that 
the elastic contribution to the ARPES spectral peaks 
near the antinodal and nodal points are determined by 
completely different scattering processes. 

Conclusions. The discovery of nanoscale inhomo- 
geneity in the local electronic structure of BSCCO- 
2212 Q 0001 has provoked an intense discussion about 
the origin of this phenomenon in cuprates and related 
correlated electron materials. In this work, we have of- 
fered strong evidence that the inhomogencity in the co- 
herence peak position is in fact driven by dopant atoms, 
located away from the Cu02 plane, whose primary effect 
on one-particle properties is to modulate the local pair in- 
teraction. This ansatz allowed us to reproduce, in model 
single-impurity and many-impurity calculations, most of 
the important correlations observed in recent STM ex- 
periments. 

The calculations reported here have been done entirely 
within a BCS framework, and as such cannot be expected 
to reproduce certain other correlations, such as the in- 
crease of the average gap with underdoping. Neverthe- 
less, we believe that our results represent an important 
step towards further modeling which may reveal the mi- 
croscopic nature of the modulated pair interaction. 
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